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The formation of inter-poly electrolyte complexes from the association of oppositely charged polymers in an electrolyte is studied. 
QnThe charged polymers are linear oppositely charged polyelectrolytes, with possibly a neutral block. This leads to complexes with 
a charged core, and a more dilute corona of dangling chains, or of loops (flower-like structure). The equilibrium aggregation 
number of the complexes (number of polycations m+ and poly anions m_) is determined by minimizing the relevant free energy 
functional, the Coulombic contribution of which is worked out within Poisson-Boltzmann theory. The complexes can be viewed 
|as colloids that are permeable to micro-ionic species, including salt. We find that the complexation process can be highly specific, 
giving rise to very localized size distribution in composition space (m+,m_). 



1 Introduction 
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2 "Electrostatic interactions are instrumental in determining the 
S .structure and function of living organisms, biopolymers and 
drug delivery systems. Charged macromolecules can self- 
assemble and aggregate into compact intermolecular com- 
plexes. This ability of oppositely charged polymers to form 
finite size complexes determines their biological function, 
which for example is important in gene transfection and com- 
^-H [pactization of DNA^~^, that provide promising alternatives to 
^ viral vectors^. Such macromolecular systems, where electro- 
^2 static forces are usually stronger than van der Waals or hy- 
\^ drogen bonds, exhibit rich behavior and structural variabil- 
V/^ ity. The structures formed by opposite charges are usually 
more stable than neutral block copolymers micelles dissoci- 
[ating upon dilution or slight change in the external conditions. 
fS| The concept of stabilization of intermolecular complexes by 
interaction of oppositely charged polymers is realized in inter- 
^ polyelectrolyte or polyion complexes (PIC) and polyion com- 
: ^ [plex micelles (PIC micelles) that can be used for drug deliv- 
ery-~-. High stability of PICs opens the possibility to use 
them as functional devices where the responsiveness to ex- 
' ternal stimuli can be connected with a function, e.g. recogni- 
tion at the molecular level ^, pH-sensitive switching devices^ 
or drug delivery carriers transporting charged objects through 
the cell membrane 

Polyion complexes have enhanced ability to undergo struc- 
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Fig. 1 Interpolyelectrolyte complexes formed by A) a linear 
polyelectrolyte (blue) and a diblock copolymer composed of an 
oppositely charged block (red) and a neutral block (black); upon 
assembling, these chains form the complex sketched, where the 
corona is made up of the neutral blocks; B) two linear oppositely 
charged polyelectrolytes (blue and red) with large asymmetry in the 
distances between charges («+ and n-). The segments with non- 
compensated charges form a charged corona of loops. 
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tural changes subject to external conditions, compared to neu- 
tral block copolymer assemblies. In addition to the response 
in change of temperature and solvent quality ^ \ the structure 
of the charged complexes can be very sensitive to changes in 
salt concentration pH^-^-^^'^^, charge ratio addition 
of ions or mixing ratio^*22. Tuning the molecular archi- 
tecture and global properties of PICs would allow for pre- 
cise control of their functional properties. Understanding the 
physics and fundamental features of the self-assembly of PICs 
is thus a challenging task. Oppositely charged polymers in 
symmetric solutions can precipitate into uniform macroscopic 
phase of polymers and ions^^'^^. The physics of aggregated 
chains of opposite charge in precipitate is somehow similar to 
polyampholites, polymers containing both positive and neg- 
ative charges dispersed along the chain—. However, if the 
distribution of charges along the chain is not random or one 
of the charged polymers is a diblock copolymer with a neu- 
tral block oppositely charged polymers can form finite size 
complexes composed of a dense polyelectrolyte core and a 
swollen corona which protects the cores from aggregation by 
steric repulsion. 

In this paper, we explore electrostatic properties and equi- 
librium structures of spherical complexes formed by oppo- 
sitely charged polymers. The stability of finite size aggregates 
results from the balance of the electrostatic attraction between 
opposite charges in the core of the complexes and the steric 
repulsion of backbone segments forming a corona around the 
core. The description of such complexes is similar to polyelec- 
trolyte micellization^^"^^, combined with thermodynamics of 
aggregation in bidisperse solution—"—. The steric repulsion 
in the corona should be strong enough to stabilize the com- 
plexes of finite size. This is possible when the hydrophilic 
blocks forming the corona of the complexes are long enough 
to oppose the electrostatic attraction in the core. In addition, 
the bare charge of the complexes can be screened by the ions 
of salt and counterions, thus changing the electrostatic forces 
and affecting the equilibrium properties of the complexes. The 
interplay of those effects will be studied considering two ge- 
ometries, as sketched in Fig. [T] 

• (case A) the building blocks are a linear uniformly 
charged polyelectrolyte and a diblock copolymer com- 
posed of a neutral block and a charged block of oppo- 
site charge. The charged blocks aggregate with the lin- 
ear polyelectrolytes to form a complex with a core sur- 
rounded by a corona of neutral segments. This case is that 
of a "hairy" and neutral surrounding outside the charged 
core. 

• (case B) the building blocks are two linear polyelec- 
trolytes of opposite charge with a large asymmetry of the 
distances between the charges. The core of such com- 
plexes is composed of charged blocks of both signs, and 



is surrounded by the corona of loops of the segments 
between the charges. The corona can be neutral (seg- 
ments between neighboring charges along the chain) or 
slightly charged (tails or longer segments between dis- 
tant charges). Compared to case A, the dangling "hair" 
are replaced by loops, that may bear an electric charge. 

The paper is organized as follows. In section |2l we first 
consider case A, with a charged core decorated by neutral 
dangling hair. The electrostatics of the complexes is taken 
into account through the full Poisson-Boltzmann (PB) equa- 
tion which is solved numerically and compared with the ana- 
lytical expression of the linearized Debye-Hiickel (DH) equa- 
tion. Particular attention will be paid to the counter-ion up- 
take, where a significant quantity of charge can be "trapped" 
inside the core, thereby reducing the electric field created out- 
side the complex. A second mechanism for charge reduction 
is ascribable to the non-linearity of Poisson-Boltzmann frame- 
work: non-linear screening effectively modifies the total core 
charge, leading in general to a reduced effective (or renormal- 
ized) quantity seen from a large distance^^. At this level of 
description, the dangling hair are not taken into account. The 
more complex situation where the charged core is surrounded 
by charged loops (case B) will be addressed in section [3] In 
turn, these results will be used in section |4]to discuss the com- 
plexation behaviour of oppositely charged polymers. Conclu- 
sions will finally be drawn in section [51 An appendix summa- 
rizes the main notations employed. 

2 Charged core surrounded by neutral corona 

(case A) 

2.1 The model and its three relevant charges 

The simplest structure of a thermodynamically stable polyion 
complexes of finite size is a spherical core, containing all bare 
charges, surrounded by a neutral corona (Figure UK). Such a 
complex can be formed, for example, by diblock copolymers 
containing neutral blocks The electrostatic in- 

teractions between oppositely charged polyelectrolytes drive 
the formation of a dense core which is stabilized by the steric 
repulsion of neutral blocks forming swollen corona around the 
core. Even in such simple geometry, it is possible to tune the 
structure of the complex. Its size can be controlled by the 
lengths of the blocks, the density of charges, pH, the charge 
asymmetry, the salt concentration and solution properties. 

Assuming that the linear polyelectrolyte is positively 
charged while the blocks of the diblock copolymer bear a 
negative charge, a linear polyelectrolyte is described by the 
number of charges on the chain z+, and the distance between 
the charges while the block copolymer is described by the 
number of charges z-, the distance between the charges ri-. 
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and the length of a neutral block, N. Hence, the length of the 
polyelectrolyte chain is and the total length of a block 
copolymer is n-Z- + N. Here, the lengths are expressed in 
units of a Kuhn length (assumed common to both cationic and 
anionic chains). 

If all polyelectrolyte charges of both signs are buried in the 
core and the neutral blocks form the corona, the total "bare" 
charge of the core formed by m+ linear polyelectrolytes and 
m_ block copolymer chains is 



Zi z+m+ - z- 



(1) 



We will assume in the subsequent analysis that this charge 
is uniformly spread over the globule of radius Rc, and there- 
fore occupies a volume 47rRl/3. The counterions and the salt 
ions can penetrate into the core of radius Rc, and thus, screen 
the bare charge of the polymers. The resulting charge of the 
"dressed" core, Z2, is the charge of the core screened by small 
ions; assuming spherical symmetry, we have 



Z2 



= 471 

Jo 



r^drp{r) 



(2) 



where the total charge density p (r) reads (in units of the ele- 
mentary charge q) 



p{r). 



47r7?3/3 



H{Rc - r) + Cooe 



-Pq(p{r) 



pq(p{r) 



(3) 



In the above relation, P = l/kT is the inverse temperature, 
<p(r) is the electrostatic potential and H{Rc — r) denotes the 
Heaviside step function, equal to 1 inside and outside the 
core. The first term on the right hand side of Eq. (|3} stems 
from the polymeric matrix, that contributes to the charge den- 
sity as a spherical uniform background. We assume here that 
the system is in osmotic equilibrium with a salt reservoir with 
equal densities Coo of labile cations and anions; the canonical 
situation, where the salt density in the system would be a pri- 
ori prescribed, is amenable to a very similar treatment as the 
one presented here. From the reservoir ionic concentration, 
we define the Debye length 1/k* through = S7i£bCoo, where 
£b = P^^/^ is the Bjerrum length and £ is the solvent dielec- 
tric permittivity. In Eq. (O, the last two terms are for the labile 
micro-ions concentration. The corresponding exponential re- 
lation between the density profiles and the local electrostatic 
potential is typical of the PB (Poisson-Boltzmann) approxi- 
mation that will be adopted in the remainder. Within such a 
mean-field simplification, the electrostatic problem at hand is 
the following: 



V2(p(r) = -f^p(r) 
^ =0 

^ =0 

dr 



(4) 



The charge of the "dressed" core Z2 should be smaller than 
the "bare" charge Zi Eq. ([T]), because of counter-ion pen- 
etration inside the globule. However, the latter charge may 
be large enough to trigger significant non-linear screening ef- 
fects, that translate into an effective (or renormalized^^'^^) 
charge Z3 that can be much smaller than Z2. To be more pre- 
cise, in the weakly coupled limit where Z2 ^ (where one 
also has Zi ^0), it is possible to solve analytically Eq. dU, 
since it reduces to V^(p = K^(p for r>Rc. The resulting DH 
(Debye-Hiickel) potential reads, for r>Rc 



(PDH{r) =Zie{KRc 



Kr 



(5) 



where is a salt-dependent geometric prefactor; the complete 
solution will be provided below in section 12.21 To define the 
renormalized charge Z3, it is sufficient to note that beyond the 
linear Debye-Hiickel regime, for an arbitrary charge Z2, Eq. 
^ again takes the form V^(p K^(p, but at large distances r 
where (p becomes small. We consequently have, within the 
non- linear PB framework: 



(p{r) ^ZsSiKRc)- 



Kr 



for 



(6) 



By construction, Z3 Z\ in the Debye-Hiickel regime, while 
Z3 <C Zi upon increasing Z\ . 

Introducing the dimensionless electrostatic potential u(r) = 
Pq(p{r) and dimensionless distance x = r/Rc, Eq. dH) in spher- 
ical polar coordinates is written in a dimensionless form 



u\0) -- 



-u(x) = -3ZiH{l -x) + {kRc)^ smh{u{x)) 



X 



= 



There are then two dimensionless governing parameters, Zi = 
Ziis/Rc and kRc. The solution of this nonlinear equation 
gives the charge density and the distribution of small ions 
around the complex. 

2.2 Two limiting cases: weak charges and salt-free situa- 
tion 

Equation (|7]) can be solved analytically in the DH approxima- 
tion when the electrostatic potential is small, <C 1. In 
this case, the charge density p (x) can be linearized, ^^"W ^ 
l±u{x), and the solution can be written in the form 



l-{l^KRc)e-^^^'-^ 
KRcCOsh{KRc) — sinh{KRc) 



r<Rc 



Kr ' 



r>Rc 

where = ?>Z\/{kRcY. As a consequence, the geometrical 
constant that enters into the electrostatic potential at large dis- 
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tances (0) is 

e{KRc) : 



[KRcCOsh{KRc) - smh{KRc)] . (9) 



Our results for a charged polyelectrolyte complex in the 
presence of salt can be compared with the salt free regime 
In this case, it is essential to enclose the complete system in 
a confining boundary, otherwise the counter-ions "evaporate" 
-their energy loss upon leaving the globule vicinity is out- 
beaten by the entropy gain of exploring a large volume- and 
the problem becomes trivial. We therefore define Rws^ the 
Wigner-Seitz radius of a large sphere containing the sys- 
tem. We note that the ratio 7] = (Rc/Rws)^ defines the volume 
fraction of globules in our system. In the present case where 
counterions only are present, the density of charges is 



(10) 



where cq is a normalization parameter to ensure total elec- 
troneutrality (it does not have any physical significance as 
such, unless a particular "gauge" or reference has been cho- 
sen for the potential). A possible choice among others is 
47rcoR^s/^ ~ -^1 • corresponding dimensionless PB equa- 
tion reads 



'-u'ix) -3Zi 

u'{Rws/Rc)=^ 



H{\-x)-rie' 



(11) 



with re-scaled distance x = r/Rc. The dimensionless charge 
Zi = Z\Ib/Rc and 7] are here independent control parameters. 

2.3 Results 

The solutions of Jhe PB equation for different dimension- 
less bare charge Zi = Zi^b/Rc are presented in Figure|2t). The 
comparison with DH approximation ([5]) shows as expected 
that this approximation is valid for weakly charged objects. 
Such a comparison is a test for the numerical procedure used 
to solve Eqs. ©, and of course, strong deviations are observed 
between DH and PB solutions for Z\ larger than a few units. 
It is noteworthy in Fig. [2] that the electric potential inside the 
core tends to a plateau when Zi is high enough. The corre- 
sponding labile ion local charge indeed tends to compensate 
for the background charge, resulting in a vanishing total local 
charge density. This requirement implies that one has 



u arcsinh[3Zi/(/c7?c)^ 



(12) 



which gives w ^ 3.91 in Fig. [2^) for Zi = 26, and likewise 
u 8.48 in Fig. [Jt) for kRc = 0.18, as can be seen. 



The solution of the salt- free equation ([TT1) is shown by a dot- 
ted line in Figure |2^). The absolute value of the correspond- 
ing electric potential u cannot be compared to its counterpart 
found with salt, but the variations of u can be. It can be seen 
on the figure, panel a), that the amplitude of u is as expected 
larger without salt (for the same value of Zi). This illustrates 
the weaker screening without salt. In addition, panel b) shows 
that the small kRc limit coincides with the salt-free limit, as 
it should (here, the salt-free solution has been shifted by the 
constant required to have the same potential at r = as in the 
kRc = 0.02 case). The salt- free results reported here depend 
very weakly only on packing fraction. 

The solution of the PB equation provides also the distri- 
bution of small ions around the core of the complex p±(r) = 
CooC^^^^^ . It is shown for different charges of the core Zi in Fig- 
ure [3^) and fixed salt concentration Coo in Figure [3]3). As the 
charge is increased the concentration of small ions of opposite 
charge inside the core increases inducing stronger screening 
effect. The concentration of ions of both signs in the core in- 
creases with increase of the bulk concentration of salt. This 
is shown in Figure [3j)), where the bare charge of the core 
Z\ is fixed, while the concentration of salt in the solution is 
changed. If the salt concentration is low, the redistribution of 
small ions around the core is almost due to counterions and 
the concentration profile of small ions obtained from Eq. (|7) 
approaches the corresponding salt free solution given by Eq. 
([TT1) . In addition, the increase of salt concentration results in 
higher concentration and induces stronger redistribution of the 
small ions around the core. 

The presence of small ions inside the complex impinges on 
the radial distribution of charges p (r). The background "core" 
contribution to this quantity is a step function, while due to 
the penetration of labile ions, p (r) is small in the globule cen- 
ter, where charge neutralization is most efficient, and increases 
upon increasing r (see Fig. |4]). In addition, the charge density 
in the center of the core vanishes at high concentrations but 
also for large enough Z\ . 

The screening of the bare charge by small ions penetrating 
into the core can be measured by the charge of the "dressed" 
core, Z2 introduced above in Eq. (O. The analytical expres- 
sion for Z2 can be obtained within the DH approximation: 
p(r) = 3Zi/(47r7?^) — 2cooUDH{^r), where UDni^r) is given 
by Eq. Thus, 



from which it follows that 
ZlDH A^kR, 



Zi {kRc^ 



{KRcCOsh{KRc)-smh{KRc))e~'^^'. 



(14) 

It can be checked that this relation is consistent with the more 
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Fig. 2 a) Dimensionless electrostatic potential w as a function of the charge of the core Zi, for a fixed salt concentration corresponding to 
kRc = 1. 11. b) Electrostatic potential w as a function of salt concentration. The charge of the core is fixed, Zi = 26. The dashed line is the 
solution of Eq. iT\\ with no salt, solid lines are the solutions of the non-linear PB equation (IT]), black dashed line is the solution of the DH 
equation ([S]), and the value of the packing fraction to solve the salt-free problem is rj = (Rc/Rws)^ = 0.000125. 



familiar DH result for the potential of a spherical colloid hav- 
ing bare charge Z2 and radius 



33. 



(p{r) = Z2 



I^kRc 



which imposes that 



^{kRc) 



I^kRc 

KRr 



(15) 



(16) 



From Eq. ([T4l) . we have that Z2 ^ Z\ within DH approxima- 
tion, up to a salt-dependent prefactor. However, upon increas- 
ing Zi , non-linear effects become prevalent and invalidate the 
DH approach, see Fig. [51 obtained by solving the non-linear 
PB theory. The corresponding slower than linear increase of 
Z2 with Zi is illustrated in Fig. [5^) and b), see also panel 
c) for the salt- free case. Increasing salt concentration screens 
out the charge, i.e. it leads to a decrease of Z2, and flattens 
the curves in panel a), where the DH prediction ([T4]) holds for 
low Zi . In essence, increasing salt concentration ultimately 
leads to the DH limit where Z2/Z1 is Z\ independent, as can 
be inferred from Eq. dS}, where w(0) is seen to decay with 
the increase of k*. A similar conclusion is drawn from expres- 
sion ([T2I) : the DH limit is reached in the high salt limit. This 
trend is clearly seen in Fig. |5j)), where all curves tend to col- 
lapse onto the DH behaviour for kRc > 10. On the other hand, 
Z2/Z1 is a strongly nonlinear function for large Z\ in the no 

salt case. More precisely, it has been shown in Ref [39| that 

1/2 

in the strongly non-linear salt- free regime, one has Z^^Z^ 



^ —1/2 

(or equivalently Z2^Z^' ). This prediction is successfully put 
to the test in Fig. [5]:), which also shows that a change in the 
packing fraction r\ = (Rc/Rws)^ does not affect the features 
discussed. 



Finally, the behavior of the charged globule at large dis- 
tances is encoded in the effective charge Z3, defined in Eq. 
(|6]), and therefore extracted from the far-field of the numerical 
solution to the non-linear equation (|7]). Such a quantity would 
rule the interactions between two distant globules. The corre- 
sponding plots of Z3 are shown in Figure |6^) as a function of 
the globule charge for fixed salt concentration, and as a func- 
tion of salt density for fixed background charge in Fig. |6j)) and 
c). As is invariably the case in such mean-field approaches, 
the effective, or renormalized, charge increases upon increas- 
ing the bare charge It also increases with salt concentra- 
tion— and as imposed by the very definition of Z3, we find that 
Z3/Z1 ^ 1 in the DH limit (enforced either from considering 
low Zi or large kRc). In addition. Fig. [St) shows that the ra- 
tio Z3 /Z2 = Z3/Z2 is independent of bare charge Zi , except at 
very small salt concentrations. This reflects the fact that even 
for large Zi , counterion uptake is such that Z2 is significantly 
reduced, and such that the colloid included internal salt ions 
can be treated by linearized mean-field theory. Indeed, it can 
be seen in Fig. ^) that Z3/Z2 is close to its DH counterpart, 
given by Z1/Z2DH, see Eq. ([T4l) . 
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Fig. 3 Reduced density (up to a factor of two) of small ions around the core p± -- 



a) for different charges of the core Zi and fixed 



salt concentration (kRc = 1. 11) b) for different salt concentrations and fixed charge of the core (Zi = 26). For the no salt case, T] = 0.000125. 




Fig. 4 Charge density pR^iB as a distance from the core for different charges of the core and fixed salt concentration kRc 
?7 = 0.000125 for the no-salt result. 



1.77. Here, 



3 Charged core surrounded by charged corona 
(case B) 

Thermodynamically stable poly electrolyte complexes can also 
be formed by the complexation of two linear poly electrolytes 
of opposite charge with a large asymmetry of the distances 
between the charges A = n-/n^ I, which can form flower- 
like structures'^. The core of such complexes with a partially 
compensated charges is surrounded by a corona of long loops 
of a polymer with a longer distance between the charges (Fig. 
[T]-B). The loops of size ri-, are neutral, but some larger loops 
and the tails can be charged. 

Thus, we can generalize the discussion of the previous 
chapter to the case where the charged core is surrounded by 
a charged corona. We assume spherical symmetry in the dis- 



tribution of the charges around the core, i.e. the charge in the 
corona depends only on the distance from the center of the 
core r. Since the charges in the loops and tails are attached 
to the core by polymer chains, the electrostatic interaction of 
those charges with the core is balanced by a weak entropic 
force due to polymer chain extension. If the electrostatic force 
is not very strong, a polymer chain carrying the charge can 
be envisioned as a Gaussian coil and the probability of radial 



3 

2na^ 



{Rc-rf 



where 



distribution of charges is P{r) ^ exp 

n is the length of the polymer chain in the corona and a is the 
Kuhn segment length. This approximation is valid for small 
charges that do not perturb significantly the statistics of the 
chains. The resulting density of charges is the sum of three 
terms: the bare charge of the core, the charge of the counteri- 
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Fig. 5 Charge of the dressed core Z2 (uptake charge) as defined in Eq. as a function of a) bare charge for different salt concentrations b) 
salt for different values of bare charge Zi c) bare charge in the salt-free case, on a log-log scale. The arrows in panel a) show the 
Debye-Huckel prediction, Eq. HM , that provides the correct limiting behaviour at small charges. Note that in panel b), the Debye-Huckel 
prediction does not depend on Zi. The no salt solution in panel a) corresponds io rf = 0.01 (which means Rc/Rws — 0.21). 




10 20 ~ 30 40 50 0,1 1 10 5 10 15 20 25 

^1 '^Fl^ 



Fig. 6 a) Renormalized charge Z3 of a globule as seen from large distances, as a function of the bare charge Zi for different salt 
concentrations, b) and c) show either Z3/Z2 or Z3 /Zi as a function of salt kRc, for different bare charges. In panel b), the continuous curve is 
for the ratio Z1/Z2 found within DH approximation, which is thus the inverse of Eq. ilAi . 



ons plus salt molecules, and the charge of the corona: 
p(r) = if(7?,-r)+c.o^-^^^W-Coc^^^^W 

+p,P(r)^-^^^^Wif (r - Rc) (17) 

where p^, positive or negative, is a parameter controlling the 
total charge of the corona and a = ± is the sign of this charge 
(+ when Pc > and — when p^ < 0). It was assumed here that 
all loops have the same length n and that ions in the corona 
are monovalent. This expression for p (r) leads to a Poisson 
equation similar to Eq. ©. 

Solution of Poisson's equation gives the radial distribution 
of the potential which is shown in Figure |7^) andlTt) for dif- 
ferent charges of the corona. The charged corona influences 
the distribution of small ions around the core, a quantity that 
is displayed in Fig. [7]:) and |7jl). Weakly charged coronas 
clearly do not modify the monotonous decrease of u with dis- 
tance that was observed in Fig. [21 but this is no longer the 



case when Pc is increased. Indeed, a point where u reaches an 
extremum [maximum in panel a) and minimum in panel b)] 
can be observed. From Gauss theorem, this coincides with the 
point where the total integrated charge over a sphere having 
the corresponding radius vanishes. The physical phenomenon 
occurring in panel a) where the charges in the corona are of 
the same sign as the bare core (assumed positive), is that the 
positive corona induces a migration of negative micro-ions in- 
side the core and its vicinity, that change the sign of the uptake 
charge Z2, which is now negative. Adding the corona charge 
to Z2, though, leads to a positive charge. Hence the charge 
inversion evidenced by the potential extremum. The density 
peak of negative micro-ions is clearly visible in panel c). On 
the other hand, when the bare core and the corona bear charges 
of opposite signs [panels b) and d)], a conjugate mechanism 
takes place: small labile cations are "sucked" inside by the 
core, which leads to an integrated charge in a running sphere 
that is positive for small spheres, and becomes negative once 
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it includes the corona. In all cases, the mechanism can be 
viewed as a corona induced local charge inversion. 

We do not repeat the full analysis of the difference be- 
tween bare, uptake, and effective charge in the present case, 
but we show how the total charge of the corona, Zcowna = 
/J r^drapcP{r)e-'^^'i'?^')H(r - Rc), depends on the parame- 
ter pc in Fig. [8l 

4 Complexation of oppositely charged poly- 
mers 

The previous sections, devoted to the electrostatics of poly- 
electrolyte complexes, have left aside the energetical aspects, 
to which we turn our attention hereafter. Once the total free 
energy of a given complex is known, it becomes possible to 
study the equilibrium behaviour, in particular the size distribu- 
tion, of an initial "soup" of individual polycations and polyan- 
ions. 

4.1 The total free energy and equilibrium complex size 
distribution 

The size of the thermodynamically stable complexes of oppo- 
sitely charged polymers is determined by the interplay of the 
steric repulsion of the chains in the corona and electrostatic 
attraction in the core. Thus, the formation of stable aggregates 
requires either long neutral blocks at least in one of the poly- 
electrolyte or a large asymmetry in the distances between the 
charges along the chain, e.g. ls. = ri- /hj^^ I. If the corona is 
composed of neutral blocks, the blocks should be long enough 
to stabilize the attraction in the core, if the corona is composed 
of loops between the charges, the segment ri- should be long 
and flexible enough to form a loop in the corona in the micelle. 

Consider then a spherical polyion complex made up of 
two polyelectrolytes of opposite charge. Each complex is 
defined by the number of polycations, m+, and the num- 
ber of polyanions, m_. If we assume dense packing of the 
monomers in the core, the radius of the core Rc can be ex- 
pressed in terms of the number of chains m+ and m_ as 

Rc = a\j^ (A^+m+ + A^_m_)] ^ , where a is the Kuhn segment, 
N± = n±q± are the lengths of the charged blocks (in units of 
a). The bare charge of the core Z\ is then also expressed in 
terms of m+ and m_, see Eq. ([T]). 

The distribution function of the polyion complexes Cm^^m^ 
is the number concentration of the aggregates with given ag- 
gregation numbers m+ and m_. The total free energy of the 
solution of polyelectrolytes of opposite charge, their counteri- 
ons and salt molecules is 



VkT 



where V is the volume of the system, /Vi+,m_ is the free en- 
ergy of the complex expressed in units of ^r, v is a molecular 
volume associated with the de Broglie length. Minimization 
of this free energy ^^'^^ with respect to Cm+,m- along with two 
conservation of mass conditions, fixing the total concentra- 
tions of polyanions, and polycations, (^+, 



P±= ^±^m+,m^ 
m±=0 



(19) 



gives the equilibrium distribution of the aggregates by their 

32 

Size 



vcm+,m- (vci,o)'"+(vco,i)'" cxp [- - m+Fi,o - m_fb,i ) ] 

(20) 

The free energy of the complex Fm+,m_ can be written as the 
sum of an electrostatic contribution, and a term accounting for 
the steric repulsion of tails/loops in the corona 

-^m+,m_ ~ ^el ~^ ^corona- (21) 

These two contributions are detailed below. 
4.2 The electrostatic contribution 

The electrostatic contribution 0.^1 is related to the semi-grand 
potential Ql^^, relevant to discuss the present situation which is 
canonical for the colloids (polymers), and grand-canonical for 
the salt entities (in osmotic equilibrium with a salt reservoir 
of density Coo). The semi-grand potential accounts for elec- 
trostatic attraction between polyelectrolytes and small ions in 
the system as well as the entropic contribution of small ions 
around the complexes. We have^S 



Pair) 



^el = J drl^^qp{r)cp{r)^kT Y^^Pair) (^In 

(22) 

where the first term is the electrostatic energy of the ionic dis- 
tribution, and the second term is the entropy associated with 
the translational movements of small ions. We note that the 
integral in Eq. (l22l) diverges for large systems (as it would 
also for neutral systems), so that we consider in the following 
the excess semi-grand potential with respect to reservoir 



m-^,m-=0 



(18) 



+ J drlc^ (23) 
Thus, the excess potential Q^/ finally takes the form 

Q„ = jdA\qp (r) (p(r) +kTY,Pa (r) In ^(24) 

J ^ Coo 

-/:r^Pa(r) + 2c. 
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Fig. 7 Electrostatic potential u{r) of the aggregate with a charged core Zi = 2, and a) a positively charged and b) a negatively charged corona. 
Distribution of small ions {KRc)^e^^ around the positively charged core, Zi = 2, surrounded by c) positively charged and d) negative charged 
corona. Here kRc = 1.77, Rc = 5a and n = 5. 



f 

Jo 



x^dx 



|ip(x)w(x) 



(25) 



For a spherical globule, this equation can be written in the 
dimensionless form as 

kT 

+p+(x) [-u{x) - 1] + p-{x) [u{x) - 1] + 2coo} 

where x = r/Rc is the rescaled distance and we used the equal- 
ity ln{p±{x)/coo) = T^W- The analysis of the free energy 
of the charged complexes suggests that the complexes with 
charged corona of the same sign as the bare core have larger 
free energy than their neutral counterparts, and thus, are less 
favorable (Figure|9t)). However, if the charges of the core and 
the corona are opposite, the electrostatic energy can be lower 
(Figure[9b)). 

The above applies for spherical globules, but leaves aside 
the particular cases (m+ = 0, m_ = 1) and conversely (m+ = 1, 



m_ = 0), where the object to be considered is no longer a com- 
plex, but a polyanion or polycation respectively. We then need 
to adapt the previous arguments to these cases of an isolated 
charged chain in a salt solution. The poly electrolyte chain 
is approximated as a cylinder of radius a with uniform lin- 
ear charge density A± oc again treated within Poisson- 
Boltzmann theory. Introducing dimensionless distance r= Kr, 
where = Stt^^Coo, the corresponding PB equation in cylin- 
drical coordinates yields, for an infinite cylinder 



1 d 
r dr 



dr 



du 
J? 
u(r 



= ± 



sinhw 

Ka 

-0 



(26) 



Here ^ = £bX± is the so-called Manning parameter (dimen- 
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Fig. 8 Total charge of the corona 'Zjcorona 
core is fixed, Zi = 2. 



/Zl as a function of PcR^iB, a) the salt concentration is fixed, kRc = 2.5, b) the bare charge of the 



sionless line charge ^^). Once this equation has been solved, 
the electrostatic contribution to and Fq^i for isolated 
chains of both signs follows from a similar calculation as that 
of Eq. (|25l): 

^ = + \ Jj'^MT) sinh(<?)) - (27) 

1 f°° 

- / rdruir) [cosh(w(?) — 1] 

2 JKa 

which was calculated per chain length N± expressed in units of 
Ib- In the following we assume that the Kuhn segment length 
of the polymer a is of the order of ^B^ Upon using the free 
energy of the infinite cylindrical macro-ion configuration, we 
neglect end effects, the consideration of which would be tech- 
nically more involved. 

Isolated chains, corresponding to m+ = 0, m_ = 1 and 
m+ = 1 , m_ = configurations are penalized by a large elec- 
trostatic energy ^ (see Figure [T0|. Indeed, these quanti- 
ties bear a large self-term, notwithstanding the solvation phe- 
nomenon, that manifests itself in the fact that ^ decrease, 
for fixed charge upon addition of salt (i.e. increase of Ka). 



4.3 Steric repulsion of loops and tails 

In the following, we consider a neutral and spherical poly elec- 
trolyte complex made up of the charged core formed by oppo- 
sitely charged poly electrolytes and surrounded by the corona. 
We consider long tails and large loops in the corona, thus, the 
corona of the complex is then approximated by a star polymer 
(Fig. [TJ-A) or a flower structure (Fig. [T}B). 



The electrostatic contribution (1281 ) is balanced by the steric 
repulsion between the tails or loops in the corona. If the 
corona consists of long neutral blocks (star polymer, case A), 
this contribution is approximated by the free energy of a star 
polymer containing m_ arms of length which is the length 
of a neutral block. This approximation is valid when the core 
is much smaller than the corona and the arms are long enough 
to use the scaling expression— 



-InA^' 



Gm_ +m_ Oi 



(28) 



In this expression, Gi are the universal exponents of the star 
polymers and their numerical values are calculated in Ref. 43. 

If the corona consist of long neutral loops and tails (flower 
structure, case B), the free energy contribution is similar to 
I), but the exponent is different. 



InA^^- 



(29) 



This exponent is calculated as follows. If the loops are formed 
by a single chain with p stickers joined together, 7c — 1 = 
G2p + 2(71 — (/?— l)Jv, where the first term is the contribution 
of the center with 2p vertices, the second term is the contri- 
bution of the two tails and the last term is the contribution of 
p—\ loops. Each loop contributes with the Flory exponent V 
in the dimension of the space d and is known numerically^^. 
If the loops are formed by m+ chains with z+ stickers and m- 
chains with z- stickers and all stickers are condensed on the 
core, the exponent is given by 

7c -1 a22^^^+2z_m_+2(Ti(m++m_)- (30) 
m+(z+ - l)dv-m-{z- - l)dv 
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Fig. 9 Electrostatic excess energy Q.gi of a PIC with a) positively charged corona and b) negatively charged corona, as a function of the charge 
of the core, Zi, for a fixed salt concentration, kRc = 7.93. 



Isolated non-aggregated polycations and polyanions are lin- 
ear polymers, thus their entropy contribution is Fcorona ^ 
-InA^^^i. Here we neglect the surface tension and hydropho- 
bic interactions between polycations and polyanions in the 
core of the complexes, thereby assuming that the electrostatic 
attraction of opposite charges is the leading contribution; hy- 
drophobic interactions may however be dominant for neutral 
complexes. 

4.4 Results 

Eq. (|2Q1) defines the equilibrium distribution function of the 
complexes c^^,m_ as a function of the geometry of the chains, 
asymmetry of the charges along the chain and salt concen- 
tration. We assume that the conformation of the polyelec- 
trolyte complex is a spherical aggregate with a core formed 
by charged blocks surrounded by neutral corona (Eigure [TJ 
A)). Since the electrostatic contribution of the core and iso- 
lated chains in the solution is the main contribution to the free 
energy, one might expect that the thermodynamically stable 
complexes would be narrowly distributed in size and have the 
minimal possible charge. Thus, the equilibrium of the free 
energy would require the compensation of the charges inside 
the core, such that the formed PIC micelles are almost neutral. 
However, in our description we allow for deviations from zero 
charge, because other contributions to the total free energy, the 
entropy of mixing, the salt concentration and the steric repul- 
sion in the corona, may shift the equilibrium. 

The distribution function of the complexes, Eq. (l2Ql) is cal- 
culated for each combination of (m+, m_) and the results 
are shown in Eig. fTTb). As an example, we plot the nor- 
malized vcyn^^ni- for a mixture of a linear polymer with the 



charge z+ = 18 and oppositely charged diblock copolymer 
with a charged block, z- = 78, and neutral block of length 
N = 200 (case A). The two polymers share the same distance 
between the charges along the chain: = ri- = l/^ = 4. 
It can be noted that all three distributions reported lie around 
the "electroneutrality line" = z-M-. In the vicinity of 

that line, the precise location of the support of the distribution 
function stems from a subtle balance of effects, as embodied 
in the free energy ([T8l ). We observe in Eig. fTTb ) that upon 
increasing the salt concentration, the equilibrium size distri- 
bution is shifted towards smaller complex sizes and becomes 
more peaked. Eor instance, for /ca = 0.3, the peak corresponds 
to m+ = 24 and m_ = 7 (for which the complex has charge 
z+m+ — z-M- = —72). A similar trend is observed while 
changing the length of the neutral block, which controls the 
repulsion in the corona. Long tails in the corona favor smaller 
complexes, and shift the equilibrium accordingly. Since the 
complexes are close to neutrality, the salt concentration mostly 
affects the electrostatic energy of free chains [Eq. (|28ll1, and 
the shift of the aggregation numbers along the electroneutral- 
ity line is mainly due to the chains in the solution. In addition, 
increasing the bulk concentration of polymers, vci^o and vco,i , 
increases the aggregation numbers, see Eig. [12] 

Eig. [13] shows the size distribution function of the com- 
plexes formed by equally charged ("matched" in terms of Ref . 
9) polymers (18,18), (44,44) and (78,78) and "unmatched" 
polymers, (18,78) and (78,18). The polymer concentrations 
are chosen in such a way that the complexes are formed close 
to the origin, which may indicate the onset of aggregation. 
Aggregation of long polymers, (78,78), occurs at smaller con- 
centrations than aggregation of short polymers, (18,18), due 
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Fig. 10 Electrostatic energy Q.^^ ^ of the rod (Eq. |28]) for the length £b as a function of a) the Manning parameter ^ (dimensionless linear 
charge) and b) salt concentration Ka. 

to the entropy of mixing, which strongly depends on the total 
length of polymers. We find that unmatched complexes [see 
the cases (18,78) and (78,18)], can also be formed if the aggre- 
gation numbers are close to the electroneutrality line (the op- 
posite charges are compensated). On the other hand, Ref ^ has 
put forward a chain recognition mechanism where matched 
cases are more prone to form large complexes, but the system 
considered there is somewhat different, involving the equilib- 
rium between three types of individual chains together with 
two and three component complexes. 

5 Discussion and conclusions 

We have developed a framework to study the formation of 
polyelectrolyte complexes from an initial arbitrary mixture of 
charged polymers, where both poly cations and poly anions are 
present in an electrolyte solution. Two situations were ad- 
dressed, as sketched in Fig. [T] for a given poly cation type, the 
poly anion is either a diblock copolymer with a long neutral tail 
(case A), or a poly anion having a different intercharge spacing 
along the backbone (case B). Coulombic attraction between 
oppositely charged polymers leads to the formation of com- 
plexes, with an a priori unknown composition. The numbers 
of chains of both types in a given complex were denoted m+ 
and m_ . These complexes were envisioned as forming hairy 
structures, where the hair/corona is either made up of dangling 
neutral chains (case A) or of loops (case B), while the core of 
much smaller spatial extension contains most of the charges 
of the polymeric backbones. We started by focusing on the 
electrostatic aspects, treated at the level of Poisson-Boltzmann 
theory. This part of the work thereby extends a previous study 
performed for salt- free systems In a second step, the re- 



sulting electrostatic free energy of the complexes was used, 
together with the entropic repulsion between tails/loops in the 
corona, to provide us with a free energy functional for an ar- 
bitrary mixture of complexes having a given size distribution 
Cm+,m- . Upon minimizing this functional under the appropri- 
ate constraints of mass conservation for both polycationic and 
polyanionic species, we obtained the equilibrium composition 
of our mixture. Whereas this optimal distribution turns out 
to give a negligible weight to configurations that depart from 
complex global charge neutrality -a property that may have 
been anticipated-, it exhibits the non- trivial feature of a high 
selectivity: out of an initial random soup of poly cations and 
polyanions, well defined complexes with precise composition 
(m+, m_) may emerge, particularly when the salt density is 
increased. 

The problem under study here is characterized by a large 
number of dimensionless parameters, and we furthermore 
made simplifying assumptions in the description, such as 
equating the Kuhn lengths for both positively and negatively 
charged polymers. We chiefly focused on the effect of chang- 
ing the salt concentration, which is an experimentally simple 
control parameter. The pH dependence of the core charge of 
the complexes has not been addressed, but it can be incorpo- 
rated for instance via the Henderson-Hasselbalch equation 
In addition, the Coulombic aspects were treated at mean-field 
level, which is adequate provided the bare charge of the com- 
plex core, Zi , is smaller than a bound of order [Rc/ {z^iB)]^^^, 
which decreases when increasing the valence z of the mobile 
micro-ions (assumed here monovalent, i.e. z = 1). Finally, we 
have neglected the structure of the core, by homogeneously 
smearing out its charge. This certainly leads to overestimate 
their free energy, due to the neglect of the corresponding neg- 



1 2 I Journal Name, 20 1 0, [vol]. 1 414) 



This journal is ©The Royal Society of Chemistry [year] 



(18,78) 




I 

I 



b). 



140 

0,1205 

0,1080 1 20 • 



100 
80 



0,05403 pp 

+ 

0,03603 gQ . 



40 
20 



(18,78) 





















































Ka=0. 


2 


































02=0.25 













10 



15 20 

m 



25 



I 
I 



30 



0,1205 
0,1080 

0,09001 

0,07201 

0,05401 

0,03601 

0,01801 

1,000E-05 



Fig. 11 Case A. a) Probability distribution function (normalized vcm+,m_) for asymmetric block copolymers of opposite charge, z+ = 18 and 
Z- = 78 with the same distance between the charges, and Manning parameter ^ =0.25. The length of neutral block (stabilizing corona) is 
N = 200, the salt concentration is Ka = 0.2, and the concentrations of individual chains is vci o = 10""* and vcq^i = 10~^. b) Projection of the 
same function on the plane (m+,m_) for three different salt concentrations. 
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Fig. 12 Projections of the probability distribution function 
(normalized vcm+,m ) on the plane (m+,m_) for the equally charged 
copolymers of opposite charges, z+ = 78 and z- = 78 as a function 
of polymer concentrations, vci^o and vco,i. Manning parameter, 
^ = 0.25 and salt concentration Ka = 0.2. The inset shows the 
points of electroneutrality of the complexes (black disks) together 
with a zoom onto the vci^o = vco,i = 10~^^ case. 




m 



Fig. 13 Projections of the probability distribution function 
(normalized vcm+^m ) on the plane (m+,m_) for different lengths of 
block copolymers of opposite charge with the same distance 
between the charges, ^ =0.25 and Ka = 0.2. Increasing 
concentrations of individual chains, vci o and vcq^i, increases the 
aggregation numbers, which move along the electroneutrality lines. 
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ative correlation energy—. 



Summary of main notations used 



q elementary charge 

a Kuhn length, assumed equal for both polycationic and poly anionic chains 

z± total charge of a chain, in units of ±^ 

n± distance between ±^ charges along a linear polymer, in units of Kuhn length 

N length of a neutral polymer block, in units of the Kuhn length 

m± number of positive/negative chains in an aggregate (core) 

Ib Bjerrum length / (ekT) defined from temperature and solvent dielectric 

permittivity 

Rc radius of a spherical aggregate/core 

Zi bare charge of a spherical core (due to polymers) 

Z2 "uptake" charge of a core (due to polymers and salt ions inside the core) 

Z3 effective (or renormalized) charge of a spherical core, relevant at large 

distances from the core center 
Z reduced charge, defined as Z = Zig /Rc 

Coo salt density in the reservoir 

Debye length, defined through = SnisCoo 
p (r) total density of charge at a distance r from core center 

Pc parameter controlling the charge of the corona (case B) 

^ Manning parameter defined as A±^b oc ig/n±, where A± is the linear charge 

of a linear polymer 
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